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*5 Abstract 
i 

Dark Matter annihilations after recombination and during the epoch of structure 
formation deposit energy in the primordial intergalactic medium, producing reion- 
ization and heating. We investigate the constraints that are imposed by the observed 
optical depth of the Universe and the measured temperature of the intergalactic gas. 
We find that the bounds are significant, and have the power to rule out large por- 

Q\ tions of the 'DM mass/cross section' parameter space. The optical depth bound is 

generally stronger and does not depend significantly on the history of structure for- 
mation. The temperature bound can be competitive in some cases for small masses 

^ or the hadronic annihilation channels (and is affected somewhat by the details of 

structure formation). We find in particular that DM particles with a large annihila- 
tion cross section into leptons and a few TeV mass, such as those needed to explain 
the PAMELA and FERMI+HESS cosmic ray excesses in terms of Dark Matter, 

. ^ are ruled out as they produce too many free electrons. We also find that low mass 

particles (< 10 GeV) tend to heat too much the gas and are therefore disfavored. 



1 Introduction 

After having settled in a state of electrical neutrality at the time of recombination, the atoms 
of the intergalactic baryonic gas had their electrons stripped off and returned to the state of 
an ionized optically thick plasma with free electrons during the process which is known as 
'reionization' [1]. The phenomenological information of how and when this happened comes 
mainly from two sets of observations, both involving the study of the absorption of photons that 
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have crossed such intergalactic medium (IGM), but of very different energies and in different 
contexts: (i) the measurement of the optical depth r encountered by photons of the Cosmic 
Microwave Background (CMB) during their journey from the surface of recombination to us 
and (ii) the analysis of the spectra of distant quasars. 

In this second class, the Gunn-Peterson test (the absence of complete absorption at Lyman-a; 
frequency) in the spectra of quasars located at redshift z < 6 allows to infer that the process of 
reionization of all the hydrogen atoms was complete by then [2] . Similarly, ultraviolet spectral 
observations suggest that helium has also been doubly ionized since redshift z ~ 3 [3]. The 
resulting free electron population between us and redshift z ~ 6 translates quantitatively into 
a value for the optical depth r of about 0.038 (see Sec. 2). 
The WM AP 5-yr measurements give, however, a higher value [4] 

r = 0.084 ±0.016 (1) 

(this number is determined combining with Supernovae and Baryon Acoustic Oscillations ob- 
servations; the CMB data alone give a very similar result: r = 0.087 ± 0.017 [5]). This implies 
that it must have existed a population of free electrons, at z > 6, that contribute to the optical 
depth the remaining 5r ~ 0.046 ± 0.016 < 0.062 (at la). 

As of now, however, the possible contributions of different kinds of high redshift astrophys- 
ical sources to reionization are still a matter of debate: whether first or second generation stars 
(Population III and II, respectively) or older quasars have given the dominant contribution, 
and the nature of its evolution with redshift is still unknown, see e.g. [6] and references therein. 

Astrophysical sources also heat up the IGM later in the evolution of the Universe: existing 
measurements of the temperature of the intergalactic medium, therefore, might also allow to 
place constraints on the kind of sources and their history, see [7, 8] and references therein. Such 
temperature is measured via observations of the Ly-a forest, caused by the small amounts of 
relic neutral gas on the emission lines of distant Quasi Stellar Objects. Observations of the Ly-a 
forest are limited to redshifts 2 < z < 5 : at higher redshifts the increasing neutral hydrogen 
density saturates the forest into a complete absorption; at lower redshifts the measurements 
become too plagued by observational uncertainties. The determination of the equation of state 
and temperature of the IGM through Ly-a observations involves also calibrations by means of 
numerical simulations, and virtually a dependence on cosmology which is taken into account 
by adding systematics to the observational uncertainties. We take as a reliable estimation the 
datapoints in Schaye et al. (2000) [8] (left panel of their Fig. 6). These points span z ~ 4.3 — > 2 
(see fig. 4b) and we will take all of them into account. To fix the ideas, we report here the 
explicit value of the highest redshift point, which will turn out to be most constraining (and is 
probably least affected by systematic uncertainties). It reads 

5.1 ■ 10 3 K < T igm < 2.0- 10 4 K (2) 

at z = 4.3 at la confidence level. 

1.1 The role of Dark Matter 

Dark Matter (DM), which constitutes about 80% of the total matter in the Universe, can 
naturally act as a source of ionizing radiation and gas heating. Indeed, most DM candidates 
are predicted to self-annihilate into Standard Model particles, thus injecting energy into the 
surrounding medium. From redshift of a few hundred down to now, DM has been gravita- 
tionally collapsing into the protohalos and halos that end up forming the cosmic structures 
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observed today; this condensation into bound systems enhances the DM self-annihilation rate 
(with respect to the contribution of the DM smooth density field) and therefore can increases 
significantly the rate of energy deposition, considering that the DM halos do not retain the 
products of the annihilation process inside the gas cocoon. 

The ionizations and heating can be produced both by the highly energetic photons directly 
emitted in the annihilation of two DM particles (with an energy of the order of the DM mass: 
typically from tens of GeV to tens of TeV) and by the lower energy photons produced by 
Inverse Compton (IC) scattering of the CMB photons on the energetic e + and e~ from DM 
annihilations. The latter turns actually out to be by far the most important process [9]: in fact, 
the cross section for scattering decreases rapidly with the energy of the impinging photon, 
so that low energy photons are more efficient in knocking off the electrons from the atoms. 1 
These 'primary reionization' electrons then deposit their energy in the intergalactic medium 
through several other interactions, freeing many more electrons (that contribute to the optical 
depth) and also augmenting the temperature of the gas. 

The effect of DM annihilations on reionization has been previously considered in a number 
of works [10, 11, 12], that mainly addressed the case of light Dark Matter candidates. More 
recently, ref. [9] has pointed out the importance for reionization of the Inverse Compton photons 
discussed above, and has studied the possibility of having DM produce the required 'missing' 
reionization 5t discussed above. 

We take the somewhat opposite approach of investigating whether the requirement of not 
exceeding the optical depth measured by WMAP can impose significant bounds on the prop- 
erties of DM. We also consider whether constraints can come from the measurements of the 
temperature of the intergalactic medium. 

The crucial parameters that determine the amount of (ionizing) energy emitted by DM 
annihilations are the mass of the DM particle m x (in first approximation, most of its rest 
mass becomes photons and heat) and the average annihilation cross section (crv). Up to one 
year ago, the benchmark values for these quantities were typically taken to be m x ~ tens to 
hundreds of GeV and (crv)thermai — 3 • 1CT 26 cm 3 /sec, the values for which the relic abundance 
of DM particles comes out, via the thermal freeze-out process, to match the observed f^M^ 2 — 
0.110 ± 0.005 [13] 2 . 

Recently, however, larger DM masses (around a few TeV up to multi-TeV) and much larger cross 
sections (of the order of 10~ 24 cm 3 /sec up to 10 _20 cm 3 /sec or more, depending on the mass of 
the candidate and the annihilation channel) have been invoked in the hope of explaining in terms 
of galactic DM annihilations 3 the anomalous signals reported in the fluxes of charged cosmic 

^or instance: for a Dark Matter particle mass of 100 GeV, the prompt photons will have a similar energy 
while the ones produced by IC scattering from CMB photons at z ~ 10 will peak at about 100 MeV; the former 
ones will have a scattering cross section on the electrons in the atoms of hydrogen of the order of 10 -29 cm 2 
while the latter ones of 10 -26 cm 2 , a difference in ionization efficiency of three orders of magnitude. We will 
discuss the relevant cross sections in detail later. 

2 Here Hdm = Pbm/Pc is defined as usual as the energy density in Dark Matter with respect to the critical 
energy density of the Universe p c — 3Hq /8ttGn, where H is the present Hubble parameter, h is its reduced 
value h = H o /W0 kms^Mpc -1 . 

3 The origin of these excesses could simply lie in ordinary (albeit possibly peculiar) astrophysical sources, such 
as one or more pulsars [14], sources of CR in galactic spiral arms [15] or secondary production and acceleration 
of CR in aged Super Nova remnants [16]. In this case, the sources would not be cosmological but instead located 
in the galactic disk, and moreover not too far from the Earth, since e ± quickly loose energy when travelling 
from farther away that about 1 kpc. These explanations will be confirmed or ruled out by further, more precise 
measurements of the spectra and possibly improved computations of the expected yields. 
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rays (CR) by several experiments: PAMELA [17], ATIC [22], PPB-BETS [23], FERMI [25] and 
HESS [26, 27], as we will discuss in detail below. 

The DM interpretation of these signals is already constrained to various degrees by other 
astrophysical and cosmological bounds: apart from the anti-proton bounds already taken into 
account [28, 29], limits come from gamma ray and radio waves from the galactic center [30], from 
ICS gamma rays from the galactic halo [31, 32], from the integrated cosmological flux of ICS 
photons [33, 34], from neutrinos from the galactic center [35], from CMB observations [36, 37] 
and from Big Bang Nucleosynthesis [38]. In particular, the works in [36] and [37] are closely 
connected with our analysis, as they also deal with the deposition of energy from DM annihi- 
lation and the impact on CMB observables. We will comment later on the complementarity of 
the studies. 

Our aim is to confront these recent DM scenarios and, more generally, the framework of 
annihilating Dark Matter (including the more traditional thermal annihilation cross section and 
tens-of-GeV mass) with the constraints that can be imposed by its role during the reionization 
epoch. We find that the bounds are indeed generally very relevant, and have the power to rule 
out large portions of the m x — (av) parameter space. 

The rest of the paper is organized as follows. In Sec. 2 we review the basic formalism which 
we use to compute the effects of DM annihilation during reionization. In Sec. 3 we discuss the 
implementation of that formalism, the choices of parameters that we adopt and the dependence 
of the results on such choices. We also present the DM scenarios that we consider and the fits 
that we perform to the CR data mentioned above. In Sec. 4 we present our results and we 
discuss the bounds that they produce. 



2 Computing the reionization and heating induced by 
DM annihilations 

In this Section we briefly discuss the formalism that we use to compute semi-analytically the 
effect of DM annihilations during the epoch of reionization. We follow in particular the discus- 
sions in Ref. [11, 12, 9]. 

The first quantity that we need to compute is the total optical depth encountered by the 
CMB photons as they travel from the surface of last scattering to us. It is given, in full 
generality, by the integral over the time of travel of the photon of the number density of 
unbound electrons n e (z) (the scattering targets) multiplied by the Thomson scattering cross 
section ctt = 8irrl/3 = 0.6652 barn (in terms of the classical electron radius r e ) 4 

n e (z)aT^-. (3) 

Here 

dt 1 1 

— = ~ (4) 

dz H (i + z)^n M (i + z)* + n A #ov / n^(i + ,2) 5/2 

is the standard ACDM relation between time and redshift, approximated for the regime of mat- 
ter domination in which we are interested. We use the central values inferred from the WMAP- 
5yr measurements [4] for the cosmological parameters today: Ho = 70.1 kmsec _1 Mpc _1 , = 
0.279, Q\y = 0.0462 and Qa — 0.721. As discussed in the Introduction, the totality of the 



4 We always work in natural units c — h = ks = 1. 
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hydrogen and helium gas is assumed to be ionized below redshift 6, and helium is also doubly 
ionized below redshift z = 3. Recalling that helium constitutes about 24% in mass [39] of the 
baryonic content of the universe (so that the number of helium atoms = 0.06 rib, while for 
hydrogen % = 0.76 r^), one can simply express r in terms of the number density of atoms 
today n A = (0.76 + 0.06) n b = 0.82 p c Q h /m p ~ 1.92 • lO^cm" 3 as 



0.88 f 3 , dt. , 3 f 6 , dt. .3' 
dz—il + zf-l dz—(l + zf 



0.82 J dz J 3 dz 



+ n A (TT 



dt 



0.038 5t 

(5) 

In the above relations, and represent the number density and energy fraction of baryons 
today (m p being the proton mass) and the factors of (1 + z) 3 rescale the densities to any redshift. 

5t denotes the amount of early optical depth caused by the unknown fraction x ion (z) of 
(singly) ionized atoms above redshift 6. Such reionized fraction obeys the differential equation 

n A (l + zf d ^l = I(z)-R(z), (6) 
or, equivalently, in terms of redshift 

-n A tf + z) 11 / 2 ^^ = I(z) - R(z). (7) 

On the right hand side are the rate of ionization per volume I(z), that tends to increase Xi on , 
and the rate per volume R(z) = Rn(z) + Rn e (z) with which hydrogen and helium atoms of the 
IGM tend to recombine even while reionization is proceeding. These recombination rates are 
explicitly given by the following expressions. For hydrogen 

76 

Rk(z) = K H n H n e - = K H:r^ ( n A(l + z) 3 x ion (z)) (8) 

where k u ~ 3.75 • 10" 13 (T igm (z)/eV) cm 3 /sec is an effective coefficient determined by fits to 
experimental data [40]. Ti gm (z) is the temperature of the IGM, also affected by DM annihila- 
tions, that we will discuss below. Similarly, for helium 

Rne(z) = K He^^ (n A (l + zfx- lon (z)) 2 (9) 

with k Hc ~ 3.925 ■ 10- 13 (T igm (z)/eV)°' 635 cm 3 /sec [40]. 

The rate of ionizations per volume produced by DM annihilations at any given redshift z is 
given by 

f mx dn 

I(z) = J dE^—(z)-P(E 1 ,z)-N ion (E 1 ) (10) 

where 4w-(z) is the spectral number density of DM-produced photons that are present at 
redshift z, which we will discuss extensively below, and one has to integrate over all photon 
energies E 1 from the H ionization energy e\ (or the He one, we here for simplicity do not 
distinguish the two) up to the maximum energy m x . P(E^,z) is the probability of primary 
ionizations per second, given by 

P(£ 7 , z) = n A {l + zf [1 - x ion (z)} ■ cx tot (£ 7 ), (11) 
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since the first terms represent the number of target atoms that can be ionized and cr to t is total 
cross section for all the interactions suffered by the DM-sourced photon and that result in the 
production of free electrons. It contains several contributions (we follow e.g. the discussion 
in [37]): the cross section for atomic photo-ionization — > e~A + [41] (dominant up to about 
1 MeV), the Klein-Nishina cross section for Compton scattering — > 7e~ [42] (dominant to 
about 1 GeV) and the cross section for pair production on matter 7 A — > e ± A' [41] (important at 
energies larger than 1 GeV). At higher energies, another processes that produces free electrons 
becomes important: pair production on CMB photons 7 7cmb —> e + e~. At redshift z < few 
hundred in which we are interested, its threshold is however above 10 TeV. We do not include 
the scatterings 7 7cmb — > 7 7, as they do not result in free electrons but just redistribute the 
photon energies. 

iVi on (i£ 7 ) is the number of final ionizations that the primary- ionization electron generated 
by a single photon of energy E y produces. It is simply given by 



-^ion (-^7) ^ion (-^ion \%) ) E^ 



nn 1 "-He 



nA e ilH n A e lHe 



Vion (-^ion (^) ) 



GeV 



(12) 



in terms of the ionization potential energies of hydrogen e^n = 13.7 eV and helium e^He = 
24.6 eV and their respective number abundances in the IGM. Here /x = 2.35 ■ 10 7 GeV -1 
corresponds to the number of ionizations that an electron of 1 GeV would end up causing if 
it were to release all of its energy in reionizations. The factor 7] lon takes into account the fact 
that only a portion of that energy actually goes into ionizations, the rest causing only heating 
and atomic excitations. Such fraction depends in turn on x\ on {z) itself, as determined by the 
detailed studies in [43, 44]: 

?7ion(a;ion(^)) = - — X g° n( ^ ■ (13) 

The spectral number density of DM-produced photons 4w-(z) present at redshift z is ob- 
tained by integrating the fluxes of photons produced at all previous redshifts (z r ) taking into 
account, with an absorption factor, the fact that some of them have already deposited their 
energies at previous redshifts. In formula? 



dn 
dK 



, , dt dN 

dz t=-(z 

x dz' dE'^ ! 



:i+z) 



l + z 



A3 



■A(z')-exp[T(z jZ \E^ 



(14) 



Here — 



z') is the spectrum of photons produced at z' by one single annihilation. The factors 

of (1 + z) 3 /(l + z') 3 rescale the number densities taking into account the expansion of the 
Universe. The absorption coefficient T reads 



t(z, z ',e: ( )^ 



dz"^n A (l + z 



ll\3 



<?x,o%{E") 



(15) 



where here E" = E' (1 + z") /(l + z'). A(z r ) represents the rate of DM annihilations per volume. 
It encodes therefore the information about the density of annihilating DM particles and in 
particular the halo formation history, that we discuss in the next subsection. 

As we already anticipated, beside producing ionization, DM annihilations have also the 
effect of heating the gas. The other important quantity that we need to compute, therefore, is 
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gm \ 



(that also enters in the recombination rates discussed above). It obeys the differential 
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equation [12] 

dT{g m (^zj 2X;g m (z) 
dz 1 + z 

1 / ^ion(^) Tqmb{z) — T igm ( Z) 2 //heat (-^ion (z))£(z)\ 

H ^hA{l + z) 5 / 2 Vl + ^ion(^) + 0.073 t c (z) 3n A (l + ^) 3 / ' 

(16) 

The first term just corresponds to the usual adiabatic cooling of the gas during the expansion 
of the Universe. It would lead to Ti gm (z) oc (1 + z) 2 . 

The second term accounts for the coupling between the IG gas and the CMB photons, that 
have a (redshift-dependent) temperature Tcmb- When the gas is hotter than the surrounding 
CMB, some of its energy is transferred to the photons and therefore the gas 'Compton-cools' 
down. On the contrary, if the gas is colder than the CMB, it is warmed up. The expression 
for the term in eq. (16) is obtained by writing the rate of change between the free electrons of 
the gas and the CMB photons as [45] dE e< _> y /dt = Aa^U k^n e {l + z) 3 (T CMB — T igm )/m e and 
then translating in terms of the rate of change of Ti gm of all particles in the gas dE e ^ — > 
3/2 &BWtot(l + z) 3 dT [gm (finally using eq. (4) to pass to redshift) [46]. In these relations 
U — £ ^~cmb i s th- e energy density in the CMB blackbody bath (with ^ the Stefan-Boltzmann 
constant [47]) and m e is the electron mass. Thus in eq. (16) t c (z) = 3m e /(8aT<;Tjk UB (z)). The 
various factors of (1 + z) 3 rescale the number densities with redshift. n e = x ion (z) Ua is the 
fraction of free electrons while n tot = n e + n u + + n H + n He = nA{xi on {z) + 1 + 0.073) contains the 
number density of all types of relevant particles in the gas, because it is assumed that collisions 
keep them at the same temperature (helium is here assumed to remain neutral, for simplicity). 

The third term accounts for the heating induced by DM annihilations. As DM injects energy 
at a rate S(z), the temperature changes at a rate given by 3/2 fc B tia(1 + z) 3 dT igra /dt = r] hea _ t S 
(then translated into a rate of change with z as usual). Analogously to eq. (13), the factor ^ heat 
expresses the fact that only a portion of the energy goes into heating. We adopt [43] 

^hcat (xionW) = C [1 - (1 - X? on ) b ] (17) 

with C = 0.9971, a = 0.2663, b = 1.3163. In terms of the quantities introduced above, the 
total energy deposited per second per volume by the photons in the intergalactic medium at a 
given redshift z reads 

r m x A n 

£{z) = / dE^ —{z) ■ n A (l + zf ■ a tot (Ej) ■ E T (18) 

Jo dtj l 

Solving numerically the coupled differential equations (7) and (16) allows to obtain two 
expressions for x\ on (z) (from which the value for St in eq.(5)) and T; gm (z), to be compared with 
the observational constraints discussed in the Introduction (eqs. (1) and (2)). We integrate the 
equations from z = 600 to z — 6. 

2.1 Structure Formation theory 

The annihilation rate per volume at any given redshift can be thought of as the sum of two parts 
A(z) = A sm (z) + A stmct (z). The former comes from a uniform density field of Dark Matter, to 
which we refer as "smooth", dominant before structure formation at redshifts z >100, and can 
be written as 

A-(z) = ^pl Mfi (l + zf, (19) 
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with m x being the mass of the DM particle, (av) the self-annihilation rate, and pdm,o is the 
"smooth" DM density today pdm,o = ^dmPc, Pc being the critical density of the Universe today. 
As DM collapses into gravitationally bound structures, the rise of local density will provide an 
increase in the rate of annihilations averaged over large volumes; such additional contribution 
from structure formation can be cast in terms of the number of halos of a given mass M to 
form at a given redshift z, and on the DM density distribution inside them, namely 



AStmCt(z) = 2^ I dM JM {z ' M) {1 + Z)3 I ^ 4 ^ 2 P-(r,M(^)). (20) 
For the halo mass distribution dn/dM we adopt the Press-Schechter formalism [48] 

dn ,* r x Pm , - ,da(R) 1 / <£? (1 + z) 2 \ . „. 

where cr(R) is the variance of the density field inside a radius R and 5 C = 1.28. 
We will consider different cases for the most common halo DM profiles Pi(r), commenting more 
about them in Section 3. The integral on the halo density squared in eq. (20) can be recast in 
terms of the virial mass of the halo 

M(z) = ^irr 3 s A vir (z)n M p(z)c 3 vil (M,z). (22) 

and the DM halo mass Mdm(#) obtained by integrating the DM profile up to the cutoff 
Cvi r (M, z) = r vir (M, z)/r s (the concentration parameter) 

O \ fC viT (M,z) 



M DM (^) = 1-^1 M(z) = 4vrr s Vs(M(z)) jf x 2 fox) dx. (23) 

Here r va is the virial radius. The integration variable is defined as x = r/r s , r s is the core 
radius of the given profile, p s (r, M(z)) = pi{M{z))j fi(x) and /«(x) is a functional form for the 
given type of profile. We discuss our choices for c v ; r (M, z) and fi(x), and their impact on the 
final results in Section 3. 

A viT (z) is the virial overdensity of the Universe due to the DM clustering at any given redshift 
(the radius within which the mean energy density in the halo is A vu .(z) times the smooth density 
at the given redshift p(z) = p c Qyi(l + z) 3 ), depends only on the given cosmology and for a flat 
ACDM universe can be written as [64] 

= n^ + 82(n M (»)-i)- 39 (n MW -i)'\ 

V Mm{z) J 

being a smooth function of the redshift. It is approximately 187r 2 for large enough redshifts. 
By defining the concentration function 



f^ (M ' z) x 2 Mxf dx 

^ fi (^) dx 



\3 

' f-c vir (M,z) 



Fi{M,z) = c vir (M,^) 3 /° x2 , (25) 



we can conveniently recast A stvuct (z) in terms of a "boost" Bi(z) due to the structure formation: 

A ■ (z) r°° dn 
Bi(z) = ^TT dMM—(z,M)F l (M,z) } (26) 
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Einasto profile 



NFW profile 



Burkert profile 




Redshift z 



Figure 1: The evolution of the effective DM density p"^ M as a function of redshift. Blue, 
magenta and orange lines refer to M m [ n =10~ 9 'Mq/ICT 6 M & /10~ 3 M Q , respectively (from top to 
bottom). The different panels assume different halo profiles. 



where M m i n is the mass of the smallest halos that form, on which we will return below. In its 
final form the annihilation rate at any given redshift reads 

A(z) = ^lp 2 DMfi (l + z) (i (l + B t (z)), (27) 

thus allowing us to define an effective, averaged DM density resulting from structure formation, 
Pim( z ) = Pdm,o (1 + <z) 3 a/1 + Bi(z) which we plot in figure 1, for different cases. We discuss it 
in the following section. 



3 Discussion 

Armed with the formalism above, we are able to compute the total optical depth and the 
final temperature of the IG gas resulting from DM annihilations. We now discuss its practical 
implementation. 

3.1 Structure formation parameters 

A critical quantity for the integration of eq. (26) is the concentration parameter c v j r (M, z), which 
can be thought of the (normalized) physical radius of a halo of given mass M. It is usually 
obtained by the results of numerical simulations, and in particular is found to be inversely 
proportional to the redshift z, namely c v ; r (M, z) =c v j r (M, 0)/(l + z) (Bullock et al. (2001) in 
[49]), as the radius of a halo of given mass grows with the redshift as the Universe expands. We 
have adopted the c v i r (M, 0) best fitting a WMAP3 cosmology [50], from [51] (Eq. 9). The core 
radius r s (M) is instead the radius of the core of a halo of given mass M, and its size depends 
on the chosen profile. In the table at page 11 we give the adopted values of r s (M) for a Milky 
Way sized halo, and the corresponding energy density p s . 

The Dark Matter profiles of the forming halos are assumed to be determined by numerical 
simulations. Recent, state-of-the-art computations seem to converge towards the so called 
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Einasto profile [60, 61] 



- (x a - 1) 
a 



a = 0.17. (21 



The Navarro-Frenck- White profile [62] and the isothermal-like Burkert profile [63] 

* wW = ft ip7? (29) ' ^W = ( i +j ki + j ») (30) 

represent instead previously standard choices. The shape of the DM profile enters in the con- 
centration function of eq. (25). We select the NFW profile as a benchmark, in which case, 
explicitly, 

FN ™(^ = f( 1 + (TT^) ( log ( 1+ ^-rT^)" 2 ' (31) 

Analogous, more cumbersome, expressions can be easily calculated for the Einasto and Burkert 
profile. One actually finds that the concentration functions for these latter profiles can be 
parameterized as F^m = ^Em-^NFW with kEm — 2.0 and Fb U t = ^Bur-^NFW with &Bur — 0.7, along 
a large range of c ww . We will present results for the three profiles. As one can already infer 
from the numerical values of the k\ coefficients, the adoption of an Einasto (/Burkert) profile 
will lead to an annihilation rate due to structures (and therefore to constraints on DM) that 
are roughly a factor of 8 stronger (/3 weaker) with respect to NFW. 

Finally, another quantity whose determination is somehow non-univocal in literature, and 
that somewhat affects the final results, is the minimal halo mass to form, M m i n , which we need 
to set the lower limit for the integration of our eq. (26). This is usually picked to be the free- 
streaming mass, Mf s , associated to a given DM particle mass, and it can be calculated from 
first principles once the type of interactions of the DM particle is assigned. An expression for 
Mf s is e.g. in [65] (Eq. 13), and the kinetic decoupling temperature can be taken from [68] 
(Eq. 19), since in the cases in which we are interested particles are typically weakly interacting 
with baryons. According to the standard paradigm, we obtain Mf s ~ lO~ 6 (/lO~ 8 /lO~ n )M for 
a WIMP mass M = 100 GeV(/lTeV/10TeV). It has been recently questioned, however, that 
the minimal halo mass is not M min = Mf s , but that the first halos collapse in fact at smaller 
redshift than previously expected, and the typical mass is M min = 10 5 — 10 7 M fs [69]. 

For our final results, we have chosen to present bounds obtained by picking M min = 1O~ 6 M ; 
we have however performed full calculations by using different values of M min , and we find that 
our constraints become stronger (/weaker) of a factor ~4 by using M min = 1O~ 9 M (/1O~ 3 M ). 

Before concluding this section, let us anticipate now that the bound from the optical depth 
will actually be insensitive to the parameters of structure formation, as it originates from the 
epoch when DM had not started collapsing in structures yet. The bound from the temperature, 
instead, will depend on the choices discussed here. 



3.2 Particle Dark Matter scenarios 

We consider, in a model independent way, Dark Matter particles that annihilate into different 
primary channels (with a 100% branching ratio). We use 

DM DM -> e + e- , /1+/J,- , t + t~ , W + W~ , 66, tt. 
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Each channel produces different spectra of prompt photons and of electrons+positrons, that we 
compute with the use of the PYTHIA MonteCarlo code [52] as described in detail in [28]. 5 We 
then calculate the spectra of IC scattering photons, produced by the energetic electrons and 
positrons from DM, as discussed in [31]. The ambient light on which ICS occurs consists here 
of the CMB radiation only, and it bears a dependance on the redshift. The sum of the prompt 
photons and the ICS photons thus computed constitutes the ^§r(z') to be plugged in eq. (14). 

We scan over a large range of DM masses m x , from light Dark Matter (10 GeV) up to very 
heavy Dark Matter (10 TeV). 

For the Dark Matter profiles of the Milky Way today p MW (r/r^ w ) we consider the three 
different models discussed above: Einasto, NFW and Burkert. The values for the parameters 
and p^ w are given by 



DM halo model 


r^ w in kpc 


pf w in GeV/cm 3 


NFW [62] 


20 


0.26 


Einasto [60, 61] 


20 


0.06 


Burkert [63] 


5 


3.15 



Note that all profiles are normalized at /?q W = 0.3 GeV/cm 3 at the location of the Earth [66] 
(see however also [67]). 

In this broad parameter space lie the regions that allow to explain in terms of DM annihi- 
lations the CR anomalies mentioned in the Introduction, as we now better detail. 

Dark Matter fits to charged Cosmic Ray data. The PAMELA satellite [17] has 
reported a significant excess of the positron fraction e + /(e + + e~) above the expected smooth 
astrophysical background and a steep rise above 10 GeV up to at least 100 GeV [18], compatibly 
with previous less certain hints from HEAT [19] and AMS-01 [20]. At the same time, no signal 
in the p fluxes has been seen, up to the maximal probed energy of about 100 GeV [21]. An 
excess in the flux of e + + e~ has also been reported by the ATIC-2 [22] and PPB-BETS [23] 
balloon experiments, that in particular found a peak at about 500-800 GeV. The ATIC-4 balloon 
flight has also confirmed its presence [24]. Later, this sharp feature has been questioned by the 
results of the FERMI satellite [25]: while an excess with respect to the expected background is 
confirmed, the e + + e~ spectrum has been found to be instead reproduced by a simple power 
law with index —3.04. The HESS Cerenkov telescope, too, has published data [26, 27] in the 
range of energy from 600 GeV up to a few TeV, showing a power law spectrum in agreement 
with the one from FERMI and eventually a steepening at energies of a few TeV. 

We perform the fits to these data 6 with the use of DM generated e + , e + + e~ and p spectra, 
as discussed in detail in [28]: we find the best fit values by scanning over the propagation 
parameters of charged cosmic rays and over the uncertainties on the slope and normalization of 
the astrophysical electron, positron and antiproton background. We do not include any galactic 
boost factor due to substructures within the Milky Way halo. 

The resulting allowed regions on the plane 'DM mass' - 'Annihilation cross section' are 
shown in fig.s 4 and 5. In order to fit the PAMELA data alone (for positrons and antiprotons), 
the DM particle has to fall in one of two classes: (a) a DM that annihilates only into leptons 

5 We stick to the case of "direct" annihilation of DM particles into a pair of Standard Model particles, 
i.e. we do not study models in which the annihilation proceeds into some new light mediator state (see for 
instance [53, 54, 55, 56, 57, 58, 32]). These "cascade annihilation processes" lead generically to softer spectra. 

6 We consider the PAMELA data for positrons at energies larger than 10 GeV only, where the uncertainty due 
to solar modulation is not present. We include the systematic error band on the FERMI and HESS datapoints. 
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(DM DM — ► e + e~, t + t~), or (b) a DM with a mass around or above a few TeV, that can 

then annihilate into any channel (i.e. DM DM — ► W + W~ , ZZ,bb,ti, light quark pairs and the 
leptonic channels above) possibly producing anti-proton fluxes at energies above those currently 
probed by PAMELA. These regions are individuated by green and yellow bands, for 95% CL. 
and 99.999% C.L., corresponding to Ax 2 (m x , (cv)) = x 2 ( m X i ( crv )) ~ Xmm — 6 and A% 2 ~ 23, 
with 2 d.o.f.. 

The addition of the ATIC (and PPB-BETS) data had pinned down the mass of the DM 
particles univocally, at about 1 TeV [28], and therefore imposed class (a) as the only possibility 
(see [28]). This is still the case with the FERMI data replacing the balloon results: the 
absence of features in the FERMI spectrum requires the DM mass to be somewhat above 1 
TeV, but the HESS indication of a cut-off limits it below a few TeV. The regions that fit 
PAMELA+FERMI+HESS combined are represented in the figures by red and orange areas 
(for the same confidence levels as above). Notice that the smoothness of the FERMI spectrum 
forbids a reasonable fit with the DM DM — > e + e~ channel, that would produce too peaked a 
feature: the reduced (Xmin)red turns out to be well above 2 for all DM profiles, so we plot no 
allowed region. 

In all cases, a very large annihilation cross section is needed: of the order of 10~ 23 cm 3 /sec 
up to 10 _20 cm 3 /sec or more, depending on the mass of the candidate and the annihilation 
channel. Our results reproduce and agree reasonably well with those of [28, 32] (for the cases 
that overlap). 

The needed large cross sections can be justified in specific models in terms of some enhance- 
ment mechanism which is effective today but not in the Early Universe, such as a Sommerfeld 
enhancement (see [70, 71, 28], and then [54, 72, 73]) due to the exchange of force-carriers at 
relative low-velocities of the annihilating particles (which has to be much more massive than 
the carrier). Although the details of the enhancement of the cross-section are model depen- 
dent, some general features can be identified: in particular, the enhancement shows an inverse 
proportionality with the relative velocity of the two particles, and it typically saturates to a 
maximum value when that is (3 < 10~ 4 . At recombination, z ~ 1000, the typical value for 
self-annihilating DM particles (with masses in the range explored by our analysis) is (3 ~ 10~ 8 , 
and the smooth density field only cools further with the expansion of the Universe. When DM 
collapses into gravitationally bound structures it soon virializes inside the halo, therefore being 
heated up by the gravity. However, most of the contribution to the annihilation power at z > 6 
comes from haloes with masses M < 1O 6 M , for which the virial velocity is approximately 
■u v ir ~ 10 km/s; this implies that almost all the annihilation contributing to our signal take 
place in environments with values of (3 ~10~ 5 , thus guaranteeing that the Sommerfeld effect is 
saturated. 

4 Results: bounds on Dark Matter properties 

Solving numerically the coupled differential equations (7) and (16) allows to obtain two expres- 
sions for Xi on (z) (from which the value for 5t in eq.(5)) and T igm (z), to be compared with the 
observational constraints discussed in the Introduction (eqs. (1) and (2)). 

In fig. 2 we show an example of how the DM-generated quantities affect these observables. 
Fig. 2a shows the evolution of the ionized fraction function of redshift for increasing 

annihilation cross sections (bottom up). This refers to a DM DM — > t + t~ annihilation channel 
and we have chosen m x = 2 TeV. One sees that, as soon as the cross section is sizable, 
the ionized fraction is already quite high at high redshifts, so it will yield a large 8r when 
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Redshift z Cross Section [cm 3 /sec] Redshift z 



Figure 2: The effects produced by DM annihilations on the observables connected to reionization 
and heating, for the case of annihilations into t + t~ . Left panel: the ionized fraction x lon as 
function of redshift for m x = 2 TeV and for increasing annihilation cross sections (erf) = (0.4, 
1.4, 5.1) ■ 10~ 23 cm 3 /sec (bottom to top). The thick red line individuates the value of the 
cross section for which the integrated optical depth exceeds the constraints. Central panel: the 
integrated optical depth 5t contributed by DM annihilations as a function of the annihilation 
cross section (crv). We plot lines corresponding to m x = (10, 70, 170, 800, 2000, 10000) 
GeV, left to right. For a given mass, values of the cross sections larger than those where the 
horizontal line of the residual optical depth St = 0.062 is crossed are excluded. Right panel: 
The temperature T lgm of the IGM as a function of redshift z, for different annihilation cross 
sections (increasing from bottom to top, the lower blue lines correspond to (0.2, 3.1) ■ 10~ 24 
cm 3 /sec). The lowermost dotted line shows the adiabatic cooling in absence of DM annihilations. 
The data points reproduce the measurements of Schaye et al (2000) [8]. In this example, any 
annihilation cross section larger than 7.9-10~ 22 cm 3 /sec (the uppermost red solid line) would lead 
to excessive heating of the IGM and is therefore excluded. The red dashed line corresponds to 
the cross section that already exceeds the 5t bound, which is therefore much more constraining 
for this example. 



integrated according to eq. (5). The red uppermost line individuates the cross section that 
exceeds 5r = 0.062, the value that fills the gap between t(z < 6) ~ 0.038 and the upper bar 
of the measured r in eq. (1) (see the discussion in the Introduction). The fact that, very early, 
the integral saturates 0.062 just means that the dominant contribution is the one from the 
annihilation of the smooth DM density field, before structures even start forming. When they 
do, at redshift z m 50, the effect is visible in the plot in a slight increase of Xi on , that however 
does not contribute much to the integrated optical depth. For this reason, we expect the bound 
from 5t to be insensitive to the parameters that govern the history of structure formation. 

Parallely, Fig. 2b shows how the optical depth 5t varies as a function of the annihilation 
cross section (on the horizontal axis) for different choices of the DM mass m x , so that here 
it is easy to see how the bound arises: for a given mass, the values of the annihilation cross 
section for which the DM-produced St exceeds 0.062 are excluded. For small masses, relatively 
small DM annihilation cross sections will already produce too much reionization. As the mass 
increases, the density of DM particles, and therefore the rate of annihilations and the energy 
injection, decreases. Therefore larger cross sections are admitted without exceeding 5t = 0.062. 

Notice that, in order to determine the bounds, we impose that DM is the only source of 
reionization earlier than z = 6. As discussed in the Introduction, it is instead astrophysically 
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Redshift z Cross Section [cm 3 /sec] Redshift z 



Figure 3: As in figure 2, but for annihilations into e + e~ and focussing on a small DM mass. 
Left panel: the ionized fraction x^ ori as function of redshift for m x = 10 GeV ((crv) = (0.2, 
0.7, 3.0) ■ 1CT 2& cm 3 /sec, bottom to top). Central panel: the integrated optical depth St. Right 
panel: The temperature T igm of the IGM as a function of redshift z. In this example, the 
maximum annihilation cross section (2.4 ■ 10~ 26 cm 3 /sec, the uppermost red solid line) leads to 
very significant heating of the IGM. The red dashed line corresponds to the cross section that 
exceeds the St bound, which is therefore a bit less constraining for this example. 



motivated to believe that quasars and early stars have contributed (it would actually be quite 
unphysical to believe that astrophysical sources did not contribute at all to reionization), but 
the size of their contribution is unknown. Including the reionizing impact of astrophysical 
sources would only reduce the room for DM, leading to stronger constraints. We also do not 
include the contribution to St from the small amount of 'relic' free electrons left over from CMB 
formation. This contribution is not well determined: if estimated conservatively at an average 
^ion.reiic ~ few • 10~ 4 it contributes (up to z ~ 700) just a small amount 5r rclic « 0.01; but, if 
larger, it could even provide a sizable portion of the required St [74], reducing considerably 
the room for DM and leading again to stronger constraints. Ours are therefore conservative 
choices. 

In fig. 2c we show how the temperature T; gm of the intergalactic medium is affected by DM 
annihilations, still for the choice of primary annihilation channel DM DM — ► r + r~, for the given 
DM mass m x = 2 TeV and for an Einasto profile. We show Ti gm (,z) for increasing values of 
the annihilation cross section (av) . At high redshifts the temperature already departs from the 
standard adiabatic cooling behavior, as the DM annihilations of the smooth component inject 
heat in the gas. At a given moment in redshift (corresponding roughly to the formation of the 
first halos), the temperature rises more abruptly. The larger the cross section, the earlier the 
effect on the temperature appears and the higher the low-redshift temperature ends up being. 
We impose to T igm not to exceed the la error bar of any of the data points in [8], reproduced 
in our fig. 4c. This individuates a maximum cross section and imposes therefore the bound. 
We thus see that this bound is regulated by the behavior of the temperature at low redshifts, 
when DM structures are forming. It will therefore depend on the choices of parameters of 
structure formation history, and for instance it would be different for a different choice of DM 
profile, rescaling with the factors discussed in sec. 3.1. For the example in the figure, that 
value of the cross section is actually already excluded by the optical depth one discussed above, 
corresponding to the dashed red line in the plot: the bound from the temperature is definitely 
subdominant in this case. 
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Given the (gently) decreasing shape of T igm at low redshifts, the highest redshift point in 
eq. (2) is often the most stringent. Analogously to the discussion above, we stress that this 
choice is a conservative one: in principle, constraints could be made tighter by choosing a lower 
allowed maximum value, as it would be the case if one assumes that astrophysical sources are 
mainly responsible for gas heating and therefore a reduced room is left for DM. 

Fig.3 shows the same quantities for the case of DM DM — > e + e~ annihilations and focussing 
on a small DM mass of 10 GeV. In this case, the ionized fraction Xi on behaves similarly to the 
case in fig. 2: it shows a more abrupt increase in the epoch of structure formation but the St is 
saturated also in this case by the high redshift part. The temperature, however, reaches higher 
values than before. Actually, the cross section for which the data points are exceeded is slightly 
smaller than the one for which St is saturated, so in this extreme case the temperature bound 
is slightly stronger. 

In fig.s 4 and 5 we show the regions on the plane 'DM mass m x ' vs 'Annihilation cross section 
{crvY that are excluded on the basis of the optical depth and temperature constraints discussed 
above, for the different primary annihilation channels that we consider and for different assumed 
DM halo profiles (used consistently for the forming halos and for our galaxy today). These plots 
are produced assuming M min = 1O~ 6 M . The main results are that: 

o Large portions of the regions individuated by PAMELA and FERMI+HESS are ruled 
out. For instance, for the case of annihilations into t + t~ or the entire PAMELA 

and FERMI+HESS region at about 2 to 3 TeV is excluded by the optical depth bound. 

o For small DM masses, the value of the thermal annihilation cross section (crv)thermai — 
3 • 10~ 26 cm 3 /sec starts to be touched by the exclusion contours. The bound from the 
temperature is most stringent for the e + e~ primary annihilation channel for the Einasto 
profile, suggesting m x > 10 GeV. For the other profiles or channels both bounds constrains 
smaller masses. 

We also comment on the following features. 

- Changing the DM profile does not affect the bounds from St, consistently with the fact, 
discussed above, that it is determined by the ionized fraction before structures even form. 
The constraint from the temperature instead scales up (/down) by a factor of ~8 (/3) 
moving from NFW to Einasto (/Burkert), consistently with the discussion in sec. 3. 
Similarly, changing the value for the minimal mass to 1O" 9 M (1O" 3 M ) would rescale 
down (up) the constraints from the temperature by roughly a factor of 4, as discussed in 
sec. 3. 

- Changing the primary annihilation channel changes the spectrum of the annihilation 
products and the fraction of the total energy (corresponding to the rest mass) that is 
injected in the form of primary annihilation electrons (that make ICS and transfer their 
energy to photons) and prompt photons. These are the only species that can be absorbed 
by the IGM and cause ionization, as discussed above. For instance, DM DM — > or 
the hadronic channels inject less energy into e + e~7 because neutrinos and protons carry 
away a part of the energy. Indeed, the bound rescales up for these channels. The precise 
ratios depend on the DM mass, the evolution with redshift and the shape of the spectrum 
(see also [37]). 
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Figure 4: T7ie regions on the parameter space 'DM mass' - 'Annihilation cross section' that 
are excluded by the reionization and heating bounds. The first column of panels refers to DM 
annihilations into e + e~ , the second into / u + and the third into t + t~ ; the three rows assume 
respectively an NFW, an Einasto and a Burkert profile. Each panel shows the exclusion contour 
due to exceeding the optical depth (blue short dashed line) and the exclusion contour imposed 
by excessive heating of the intergalactic gas (red long dashed line). We also report the regions 
that allow to fit the PAMELA positron data (green and yellow bands, 95 % and 99.999 % C.L. 
regions) and the PAMELA positron + FERMI and HESS data (red and orange blobs, 95 % and 
99.999 % C.L. regions). The horizontal orange band indicates the typically preferred value for 
the thermal annihilation cross section. 
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Figure 5: As in fig. 4, but for W + W~ , bb and tt annihilation channels. Since a DM particle 
fitting the PAMELA data has to be multi-TeV for these channels, the green/yellow bands are 
confined to large masses. There is no possibility to fit the FERMI and HESS data in these 
channels. The vertical cut indicates the kinematic threshold for the production of the primary 
annihilation particles. 
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- The reionization bounds from 5t scales approximately with (av)/m x , for all channels and 
profiles. This behavior is not of immediate understanding, but can be traced back to the 
energy that DM annihilations inject in the IGM oc n^ M (crv) m x = p^ M (av) / m x , to the 
dependence of the spectrum of e + e~7 from annihilation on the energy (to be rescaled by 
redshift) and so ultimately to the fact that the integral in eq. (14) bears a global weak 
dependence on the redshift during the smooth DM period. When instead structures start 
forming, the dependence on z in the boost factor modifies this scaling behavior. The 
bound from the temperature (that originates from this later phase, as said) scales more 
steeply on the plane m x — (av) . In the case of the hadronic channels and for not too large 
masses, however, the different shape of the e + e~j spectra partly compensates, and the 
scaling resembles the one linear with (av)/m x . 

5 Conclusions 

We have computed the constraints on Dark Matter annihilations imposed by the reionization 
and the heating of the intergalactic medium. 

We have calculated the flux of energy injected from DM annihilation, which results in ioniza- 
tion and heating of the intergalactic medium, including the important effect of Inverse Compton 
Scattering photons produced by the energetic e from DM on the CMB. We have followed the 
evolution of the smooth DM density and the formation history of DM halos from redshift of 
a few hundred to today, and solved the evolutions of the population of free electrons (which 
determines the optical depth of the Universe r) and of the temperature of the intergalactic gas 
T; gm , comparing them with the respective observational measurements. 

Fig.s 4 and 5 show our main results, in terms of excluded regions on the plane 'DM mass' 
vs 'Annihilation cross section'. We have considered several primary annihilation channels and 
several choices for the DM distribution profiles (consistently used for the forming halos and for 
our galaxy today). 

We have found that large portions of the regions that allow to fit the PAMELA and 
FERMI+HESS CR excesses in term of DM annihilations are ruled out by the optical depth 
bound. For instance, the entire PAMELA and FERMI+HESS region for the t + t~ or fi + fi~ 
case around m x = 2 TeV is excluded. 

We also found that DM particles with small masses, of the order O(10 GeV), tend to 
produce too much heating even with thermal annihilation cross sections ~ 3 • 1CT 26 cm 3 /sec). 7 
The bounds are most stringent for the leptonic primary annihilation channels and for an Einasto 
profile, suggesting m x > 10 GeV in case of an e + e~ channel. They become weaker for other 
annihilation channels (as less energy is injected into radiation that can be absorbed by the 
IGM) and DM profiles. 

In all the parameter space that we studied the dominant contribution to the optical depth 
is generated by the annihilations of the smooth DM field, before structures even start forming, 
and so the constraints on DM from this quantity are essentially insensitive to the uncertainties 
on the structure formation parameters. The increase in the temperature originates instead 
mainly from annihilations during the period of DM clustering, so the constraints are somewhat 
sensitive to the history of such clustering. The bound from the optical depth is generally 

7 This is in analogy with the findings of previous works [10], where however DM masses of the order of MeV 
were under consideration. The inclusion of the ICS photons from DM increases here the energy deposition by 
several orders of magnitude and so increases significantly the DM mass at which there can be an impact. 
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stronger but the one from the temperature can become competitive for small masses and for 
the hadronic channels, in the case of an Einasto profile. 

We stress that these constraints are produced by imposing to the 5r and to the T igm induced 
by DM not to exceed the la error bars of the respective measurements and are derived assum- 
ing no competing contribution of optical depth and heat from ordinary astrophysical sources. 
Including these contributions, with a specific prescription for their impact, would reduce the 
room for the DM contribution and make the bounds much stronger. 

The results of our analysis are compatible with Galli et al. [36] who included the contribution 
of the only smooth DM density field; their MonteCarlo analysis on several CMB parameters, 
clearly dominated by the constrain on r from the WMAP5 data, is however parametrized in 
terms of an efficient coupling between the DM shower induced and the IGM gas, /. In such 
formalism, the parameter / has been computed for a wide range of secondaries by Slatyer et 
al [37], that included all the relevant processes we have self consistently implemented in our 
analysis. All the results are consistent. The actual bounds agree within a numerical factor of a 
few. Here we have presented a separate estimate of the contribution to r given by the smooth 
and structure boosted annihilations, finding in particular that the smooth one is dominant. 

For the regions that are still below the bounds but where the effect of DM is already sizable, 
it would be interesting to investigate whether the resulting ionized and warm intergalactic gas 
can modify the process of the formation of stellar objects. This kind of analysis is however 
beyond the scope of our study. 

Finally, notice that we have not treated the case of decaying Dark Matter, which has 
been also recently proposed as an interpretation of the cosmic ray data (see [75] for model 
independent studies and for references). In this case, the relevant quantities are the DM mass 
and lifetime (as opposed to the annihilation cross section) and the DM energy injection is 
proportional to the first power of the DM density /?dm instead of the square. Effects are 
therefore expected to be generically milder. Ref. [76], which studies the contribution of the 
smooth DM density field, finds lower bounds on the lifetime that remain somewhat below the 
values required by current CR data. For a diffuse contribution like the one we are studying, and 
in the case of decaying DM, the energy injection averaged over the whole Universe depends on 
the total mass content, and will therefore be unaffected by clustering due to structure formation. 
A dedicated study would however be necessary to address this issue more quantitatively in our 
formalism which, as opposed to the one in Ref. [76], treats explicitly the coupling between the 
DM-induced shower and the IGM gas. 

Note Added. As this work was in the final stages of completion, on the arXiv appeared 
ref. [77], that also discusses the bound from the optical depth, for leptonic DM annihilation 
channels and masses above 100 GeV. Our conclusions are consistent and in agreement with 
theirs, for the cases that overlap. 
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